Modeling spin transport with current-sensing spin detectors 
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By incorporating the proper boundary conditions, we analytically derive the impulse response (or 
"Green's function") of a current-sensing spin detector. We also compare this result to a Monte-Carlo 
simulation (which automatically takes the proper boundary condition into account) and an empirical 
spin transit time distribution obtained from experimental spin precession measurements. In the 
strong drift-dominated transport regime, this spin current impulse response can be approximated 
by multiplying the spin density impulse response by the average drift velocity. However, in weak 
drift fields, large modeling errors up to a factor of 3 in most-probable spin transit time can be 
incurred unless the full spin current Green's function is used. 



Measuring spin-polarized electron transport requires 
both successful spin injection and spin detection. [l|, H HJ 
To extract meaningful spin transport parameters, mod- 
eling the response characteristics of these components is 
essential. In particular, the details of the spin detection 
mechanism can impose very different boundary condi- 
tions on the spin distribution and affect the device output 
response. 

Two types of electrical spin-detection techniques are 
currently in use: voltage sensi ng 3, M , d, 0, Hi and cur- 
rent sensing^ El El El El 1 14L Il5j. Because voltage 
sensing (which is sensitive to the spin density) uses an 
open-circuit configuration, ideally it does not serve as a 
sink for spins and therefore the presence of the detector 
does not impose any boundary conditions on the spin 
distribution. However, current sensing methods (which 
are sensitive to spin current flowing into the contact) are 
ideally perfect sinks of charge and spin, and so do indeed 
impose a boundary condition. Specifically, they should 
constrain the spin density to vanish at the detector con- 
tact, so the functional response of the two detectors will 
be very different. Here, we show how this boundary con- 
dition can be incorporated into the problem to accurately 
model current-sensing spin detectors, and find the regime 
under which large simulation errors are incurred by ig- 
noring it. 

To model the spin detector, the transport impulse re- 
sponse (or "Green's function" , i.e. spatio-temporal evo- 
lution of spin density with initial conditions s(x, t = 0) = 
5{x), where S(x) is the Dirac-delta function) is required. 
The spin density s(x,t) for < x < L (where L is spin 
transit length; the injector is at x = and the detector at 
x = L) and t > is determined by the spin drift-diffusion 
equation 
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function, we can construct the response to arbitrary spin 
injection conditions, even the DC currents used so far in 
experiments. 

With voltage-sensing spin detection, there are no 
boundary conditions on s(x, i) (besides those at infin- 
ity which keep the solution bounded). Assuming v is a 
constant independent of position x, the Green's function 
of Eqn. [T]can be determined straightforwardly using e.g. 
separation of variables and the Fourier representation of 
the Dirac delta function. It is given by 
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The impulse response at the detector can be obtained 
simply by substituting x = L. An example spin density 
distribution using Eqn. [2] and v — 10 6 cm/s, L — 10/xm, 
and D = 200cm 2 /s is shown in Fig. [TJ where the effects 
of spin relaxation with finite r are ignored. 

However, the Green's function of Eqn. Q] in Eqn. 
is not consistent with current-sensing detection. Un- 
like in open-circuit voltage spin detection, electrons cross 
into the detector contact. Once they cross the transport 
channel boundary, they do not return via diffusion, espe- 
cially in ballistic hot-electron spin detectors where inelas- 
tic scattering prevents it@, El El 03, El El El- This 
imposes an absorbing boundary condition on the spin 
density at the detector (s(x = L,t) = 0), from which 
spin current is determined by Fick's law, 



ds 

J s (x =L,t) = -D — \ X=L . 

ax 



(3) 



where v is the drift velocity, D is the diffusion coefficient, 
and r is the spin lifetime. Once we calculate this Green's 



The Green's function satisfying this boundary condi- 
tion can be found by making an ansatz equivalent to the 
method of images so that we can construct the solution 
using Eqn. [21 The negative-valued "image" spin density 
also moves at drift velocity v but its relative magnitude 
is determined by the transport parameters such that the 
boundary condition s(x — L, t) — is satisfied: 
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FIG. 1: Comparison of the time dependence of spin density at 
the detector using a voltage-sensing contact, spin current us- 
ing a current sensing contact (which imposes a boundary con- 
dition on spin density), and a direct Monte-Carlo distribution 
which automatically imposes the same boundary condition. 
Parameters v — 10 6 cm/s, D — 200cm 2 /s and L = 10/mi have 
been used to emphasize the relative differences, which are 
proportional to 2D / Lv to first-order and therefore decrease 
rapidly for higher drift velocity and longer transit length. 



FIG. 2: Evolution of the Green's function of the spin drift- 
diffusion equation, where a Dirac delta spin density distri- 
bution is injected at x = at t = 0, for time increments 
shown in the legend. Drift velocity is 10 6 cm/s and diffusion 
coefficient is 200cm 2 /s. For clarity, no spin depolarization 
over this timescale is considered. Note the effect of absorbing 
boundary conditions at the detector (x = 10/mi). The corre- 
sponding spin current distribution, given by Eqn. O is shown 
in Fig. [1] 



The corresponding spin current at the detector, obtained 
using Eqn. [3l is then 
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The time-evolution of Eqn. [4] (using the same values for 
v, D, and L as in Fig. [JJ) is shown in Fig. [5] The 
corresponding spin current, obtained by using Eqn. [3] 
and given by Eqn. [SJ is compared to the spin density 
from Eqn. [2] in Fig. [TJ Significant differences between 
the qualitative form and position of the two functions are 
evident with these particular transport parameters. 

Note that the expression in Eqn. [5] is simply the spin 
density in the absence of the boundary condition (Eqn. 
[2]), multiplied by the spin velocity L/t. This is a particu- 
larly satisfying result in light of the normal definition of 
particle current. 

Because current-sensing measurements are typically 
under conditions of str ong drift fields (where v >> 
D/L)[! El El El El H El, we have approximated 
this L/t term with a constant (e.g. the average drift 
velocity v) in previous work comparing experiment to 
theory. [H EH, El This approximation is justified in the 
strongly drift-dominated regime, because it involves only 
a rigid shift of the most probable transit time from 
w L/v-D/v 2 + 2D 2 /v 3 Lto « L/v -3D/v 2 + 9D 2 /2v 3 L, 
a difference of 2D jv 2 to first order and a subsequent er- 
ror in the measured velocity of 2D/L. Because typical 
values for these variables are on the order of v > 10 6 
cm/s, L > 10~ 2 cm, and 10 2 < D < 10 3 cm 2 /s, the rela- 
tive error associated with this approximation (2D /Lv in 
both cases) is then just a few percent. 



However, the error can be substantial at low drift ve- 
locity when diffusion is the dominant transport mech- 
anism. In the case that D >> Lv, use of the den- 
sity Green's function (Eqn. [2]), without incorporation 
of the functional dependence of the L/t term, will give a 
calculated most-probable (peak) transit time of L 2 /2D, 
whereas the correct spin current distribution (Eqn. [5]) 
will give a value three times smaller. This large potential 
error highlights the importance of using the appropriate 
Green's function to model the spin transport, even when 
weak drift fields are used. In addition, it should be noted 
here that in diffusion-dominant transport, the current- 
sensing detector sees a response which has a smaller stan- 
dard deviation At, resulting in smaller spin dephasing 
and therefore higher spin coherence than voltage-sensing 
detectors. 

A consistency check of this spin-current impulse re- 
sponse derived above can be provided by direct model- 
ing of the drift and diffusion of an ensemble of electrons. 
This simulation can be used to assemble a histogram of 
transit times, which in the limit of large numbers of elec- 
trons yields the transit-time distribution function. Such 
a "Monte Carlo" method, which incorporates the proper 
boundary conditions for our current-based spin detection 
technique automatically (because electron transport sim- 
ulation ceases once x > L), was already used to model 
spin transport through doped Si [1 71] . where Eqn. [TJ no 
longer has constant cocfficcnts due to an ir-depcndent 
drift velocity caused by inhomogeneous electric fields 
from ionized dopant impurities. 

The Monte-Carlo technique involves discretely step- 
ping through time a duration St, modeling drift with a 
spatial translation vSt, and diffusion with a translation of 
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FIG. 3: Comparison of most-probable (peak) and average 
(mean) spin transit times calculated using spin density with- 
out boundary conditions (Eqn. [2J, spin current with absorb- 
ing boundary conditions (Eqn. [SJ, and using a Monte-Carlo 
method. Transport parameters are the same as those used in 
Figs. 1 and 2. 

±V2D6t (where the sign is randomly chosen to model the 
stochastic nature of the process), until x > L. Then, the 
total transit time is recorded. This process is repeated 
many times until the list of transit times can be used to 
make a reliable histogram. 

A simulated Monte-Carlo distribution using 4 x 10 4 
electron drift-diffusion trajectories and St = 10 s is 
compared to the analytic functions of spin density and 
spin current (using the same transport parameters v, D, 
and L) in Fig. [TJ Clearly, the Monte-Carlo distribution 
closely resembles the spin current distribution, confirm- 
ing that this method does indeed automatically incorpo- 
rate the appropriate absorbing boundary conditions. 

Fig. [3] shows the drift velocity dependences of the 
most probable (peak) and average (mean) transit times 
of the spin density, spin current, and Monte-Carlo distri- 
butions using the same transport parameters as in Fig. [T] 
Clearly, the differences between the correct (spin current 
and Monte-Carlo) impulse response distributions and the 
spin-density distribution becomes rapidly smaller as v in- 
creases. It should also be noted that for transit lengths L 
larger than the lO^tm used here, the differences are like- 
wise smaller due to the functional form of relative error 
(D/Lv). 

As has been shown elsewhere, the real part of 
the Fourier transform of spin precession data mea- 
sured in a perpendicular magnetic field reveals the em- 
pirical transit time distribution. [181] We can therefore 
compare experimentally-obtained distributions to the 
analytically-derived models here with data from 10/im- 
thick undoped silicon spin transport devices utilizing 



current-sensing spin detectors [l,[Tl|, as in Fig. [4] Exper- 
imental results shown in the inset are obtained with an 
internal electric field of approximately 50V/cm at a tem- 
perature of 90K. Because of a rapidly degrading signal-to- 
noise, this is nearly the lowest electric field and hence the 
weakest drift conditions that can be reliably measured 
with these devices. Fitting parameters used to produce 
the analytic spin current (Eqn. [S]) and spin density (Eqn. 
[2]) are v = 5.5 x 10 6 cm/s and D =120 cm 2 /s (consistent 
with D for measurements at higher fields). The empir- 
ical distribution closely matches the spin current distri- 
bution, again confirming the importance of using Eqn. [S] 
in modeling current-sensing spin detectors. 

In summary, we have derived the correct drift-diffusion 
Green's function for modeling the response of current- 
sensing spin detectors and favorably compared it to 
Monte-Carlo and empirical spin transit time distribu- 
tions. In the strongly drift-dominated regime, the cor- 
rect Green's function can be approximated by ignoring 
the boundary condition at the spin detector and simply 
multiplying the spin density Green's function by the av- 
erage drift velocity v. However, in weak drift fields or in 
conditions of strong diffusion, large modeling errors can 
be made by using the wrong distribution. 

We thank B. van Wees and T. Stanescu for important 
discussions on this topic. 
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FIG. 4: Comparison of empirical spin current transit time 
distribution with spin density calculated without boundary 
conditions (Eqn. [2J) and spin current with absorbing bound- 
ary conditions (Eqn. [5j|. Inset: Spin precession data used to 
calculate empirical distribution with Fourier transform. [3] 
Shown is the average of four perpendicular magnetic field 
sweeps, symmetrized to avoid ferromagnetic contact magne- 
tization switching and constrain the transform to be real. 
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